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2 A minimalist stochastic model of primary soil salinity is proposed, in which 

3 the rate of soil salinization is determined by the balance between dry and 

4 wet salt deposition and the intermittent leaching events caused by rainfall 

5 events. The long term probability density functions of salt mass and concen- 
e tration are found by reducing the coupled soil moisture and salt mass bal- 

7 ance equation to a single stochastic differential equation driven by multiplica- 

s tive Poisson noise. The novel analytical solutions provide insight on the in- 

9 terplay of the main soil, plant and climate parameters responsible for long- 

10 term soil salinization. In particular, they show the existence of two distinct 

11 regimes, one where the mean salt mass remains nearly constant (or decreases) 

12 with increasing rainfall frequency, and another where mean salt content in- 

13 creases markedly with increasing rainfall frequency. As a result, relatively 

« small reductions of rainfall in drier climates may entail dramatic shifts in long- 
is term soil salinization trends, with significant consequences e.g. for climate 
i6 change impacts on rain-fed agriculture. 
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1. Introduction 

17 Large areas of cultivated land worldwide are affected by soil salinity. Szabolcs [1989] 

is estimates that 10% of arable land in over 100 countries, and nine million km 2 are salt 

w affected, especially in arid and semi-arid regions [Tanji, 1989]. Salinity refers to large 

20 concentrations of easily soluble salts present in water and soil on a unit volume or weight 

21 basis (typically expressed as electrical conductivity (EC) of the soil moisture in dS/m, 

22 i.e. deciSiemens per meter at 25° C; for NaCl 1 mg/1 ~ 15 ■ 10~ 4 dS/m). High salinity 

23 causes both ion specific and osmotic stress effects, with important consequences for plant 

24 production and quality. Normally, yields of most crops are not significantly affected if 

25 EC ranges from to 2 dS/m, while above levels of 8 dS/m most crops show severe yield 

26 reductions [Ayars et al, 1993; HUM, 2000]. Prevention or remediation of soil salinity 

27 is usually done by leaching salts, and has resulted in the concept of leaching require- 

28 ment [Richards, 1954; HUM, 1998; Schleiff, 2008]. Alternative amelioration strategies by 

29 harvesting salt-accumulating plants appear to be less effective [Qadir et al, 2000]. 

30 Salt accumulation in the root zone may be due to natural factors (primary salinization) 
si or due to irrigation (secondary salinization). Several detailed numerical models have been 

32 developed to model soil salinization [e.g. Eldin et al, 1987; Schoups et al, 2006; Corwin 

33 et al, 2007]. Generally, these models simulate unsaturated soil water flow via Richards 

34 and solute transport equations. These models are more suitable for local and short-term 

35 simulations, as they require precise soil characterization and are computationally demand- 

36 ing. Moreover, it is often difficult to identify cause-effect relationships or to synthetically 

37 compare the effects of different parameter scenarios from their numerical simulations. 
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as Vertically-averaged soil moisture and salt balance equations have also been used [Allison 

39 et al, 1994; HUM, 2000]. Despite their simplicity, these models have the advantage of 

40 parsimony, thus allowing a direct analysis of the interplay of the main processes, and 

41 provide an ideal starting point to include external, random hydroclimatic fluctuations 

42 in the analysis of long-term salinization trends. The goal of this Letter is to offer a 

43 first step in this direction. With this purpose, here we present a minimalist model of 

44 soil primary salinization, describing analytically the long-term dynamics of salt in soils 

45 caused by wet (rain) and dry (aerosol) deposition. Our aim is to quantify the salt mass and 

46 concentration probability density functions (pdfs) in the root zone, and the probability 

47 of crossing the crops salt tolerance threshold as a function of the main hydro-climatic 

48 parameters. The model framework is potentially extendible to systems including salt 

49 input from groundwater and irrigation. 

2. Methods 

50 Our starting point is a spatially lumped model [Bras and Seo, 1987] for the vertically 

51 averaged dynamics of soil moisture and salt in the root zone. As a first step we will not 

52 consider input of salt due by irrigation or groundwater upflow. Following Rodriguez- Iturbe 

53 et al. [1999], Rodriguez- Iturbe and Porporato [2004] and Porporato et al. [2004], rainfall 

54 (R t ) is modeled as a marked point process with frequency A p and with daily rainfall 

55 depths exponentially distributed with mean l/j p . The averaged soil moisture dynamics 

56 are modeled assuming constant (spatially and temporally averaged) soil and ecohydrolog- 

57 ical parameters, i.e., root depth, Z r , porosity, n, and maximum evapotranspiration rate, 
sa ET max . Assuming a rain salt concentration Cr and a constant input of salt mass per 
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59 unit ground area and per unit time by dry deposition, the root-zone mass balance for soil 

60 moisture and salt mass m is given by: 



61 where C is the salt concentration in the root zone; L(s) represents deep percolation, while 

62 ET(s) represents the losses resulting from plant transpiration and soil evaporation. As 

63 in Porporato et at. [2004], ET(s) is assumed to be linear in the range of soil moisture 

64 comprised between the wilting point, s w , and a suitable soil moisture threshold Si (an 

65 effective field-capacity threshold), at which ET occurs at the maximum rate ET max . All 

66 the rainfall input that cannot be accommodated is assumed to be lost as L(s) at si. In 
e? this minimalist model the effect of salt-induced changes in osmotic potential may only 

68 indirectly be taken into account through an average reduction of ET max . This is simply 

69 done here by keeping the same ET max for the two models (previous studies [Viola et al, 

70 2008] have shown that, in the absence of osmotic effects, the minimalist model should have 
?! artificially higher ET max to account for percolation losses below si). A full account of how 

72 reduction in evapotranspiration affects salinization patterns (reduced evapotranspiration 

73 in turn increases the available soil moisture and thus reduces the concentration of salt in 

74 the soil and increases leaching frequencies) will be given elsewhere. 

75 A complete numerical model, in which the impact of osmotic stress in reducing ET is 

76 explicitly included [Bras and Seo, 1987], has been also studied. Moreover, in the detailed 

77 model runoff takes place at saturation (s = 1), while percolation occurs for s > s/ c (the soil 

78 moisture field capacity), and it is proportional to the soil hydraulic conductivity K sat s c , 



nZ r — = -ET(s) - L(s) + R t , 
^ = C R R t +M d -CL(s), 



(1) 



(2) 
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79 where c is a soil-pore connectivity index and K sat is the saturated hydraulic conductivity 

so [Rodriguez- Iturbe and Porporato, 2004]. A comparison between the results of the two soil 

si moisture models, presented in Figure la, suggests the viability of the simplified model. 

82 Simulations for wetter climates confirm this result. 

83 The system (1) and (2) can be further simplified if one considers that the typical 

84 timescales for salt mass dynamics in the root zone are orders of magnitude larger than the 

85 ones characterizing rainfall (and thus wet deposition). Moreover, soil moisture typically 
se reaches steady-state conditions within a growing season (e.g., < 5 — 7 months), while the 
87 salt mass balance only does so on much longer times scales (e.g., > decades). Accordingly, 
ss at those long timescales, say T, the salt mass input flux can be assumed to take place 
m at a constant rate, T, that is J t t+T (M d + C R R t )dt' ~ M d T + TC R \ p /-f p = TT, and be 
oo interrupted by instantaneous and unfrequent leaching events induced by percolation. As 
9i a result, (2) can be rewritten as 



92 Leakage may then be modelled as a marked point process, with percolation depths ex- 

93 ponentially distributed with parameter 7 P [Botter et al, 2007]. For reasons of analytical 

94 tractability, the percolation events are assumed to occur according to a Poisson process 

95 with frequency A given by the frequency of soil moisture crossing the threshold s = s±. 

96 This can be expressed in terms of the soil moisture pdf as A = p(si)p(si), where the 

97 term p(s) = (ET(s) + L(s))/nZ r represents the normalized catchment- scale loss function 

98 (i.e. the total losses from the system due to evapotranspiration and leakage as a function 

99 of the soil moisture) [Rodriguez- Iturbe and Porporato, 2004]. Adopting the soil moisture 




(3) 
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100 minimalist model (for which the pdf is a truncated gamma distribution, e.g., Porporato 

101 et al. [2004]), the leaching frequency is A = n exp(— ^y^p^ /T(\ p /rj, 7) [Botter et al, 

102 2007], where T(x, y) is the lower incomplete gamma function, n = ET max /(nZ r (si — s w )) 

103 and 7 = 7 P nZ r (si — s w ). 

104 A leaching-efficiency parameter b is used to account for incomplete salt dissolution, 

105 further assuming that the typical value of soil moisture during leaching events can be 
wo approximated by the value si. With the above assumptions, the dynamics of the salt 

107 mass in the root zone can be described by a single equation 

£ = t-^, (4, 

108 where L' t is a marked Poisson noise [Van Den Broeck, 1983] with frequency A, and (di- 
ms mensionless) exponential marks with mean 

b 

n Z r si 7 P 

Figures lb and lc compare the results of both salinity models. The free parameters 
Si and b are fitted with respect to the complete model of salt mass and concentration, 
respectively. 

3 From a mathematical viewpoint, equation (4) is a stochastic differential equation with 
multiplicative white (jump) noise. In our case, since the soil solution can be considered in 
5 equilibrium during leaching events, one has to interpret (4) in the Stratonovich sense [ Van 
e Den Broeck, 1983]. Accordingly, the normal rules of calculus are preserved, and equation 
(4) can be transformed into 
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us where y(t) = \n[m(t)}. 

3. Results and Discussion 

iw The stationary solution of (6) can be obtain as in Rodriguez- Iturbe et al. [1999]. Then 

120 using the derived distribution for m, i.e., p(m) = p(y)dy /dm, we obtain the probability 

121 distribution for the salt mass in the root zone 

p(m) = A/" exp(-^r) m 1/M , (7) 

122 where J\f = (A/T) - ^/!^^^) and T(x) is the Gamma function. Equation (7) summarizes 

123 the soil salinity statistics as a function of climate, soil and vegetation parameters. 

124 Figure 2 is a graphical representation of the dependence of the mean salt concentration 

125 (C) = (m)/nZ r {s)oia the yearly rainfall and A p . The contour-lines connect equal values 

126 of the mean salt concentration in the soil, for a given input of salt T. The latter has been 

127 calculated for two different geographic regions. Typical salt inputs in coastal areas are 

128 100 — 200 kg/(ha yr) of salt, while values drop of an order of magnitude in continental 

12 9 regions [HUM, 2000]. 

uo Between the black region and the light gray ones in Figure 2a, the behavior of (C) 

131 changes substantially. Above a certain total rainfall per year, the input of salt related 

132 to rainfall frequency becomes immaterial as leaching effectively washes out the salt mass 

133 from the root zone. For lower total rainfall values, however, the salt in the soil increases 

134 with increasing A p . For a given annual precipitation depth, with low rainfall frequencies, 

135 rainfall events carry enough water to trigger leaching. Conversely, if A p is high, evapo- 

136 transpiration dominates, leaching is largely reduced, thereby causing salt accumulation 
is? in the root zone. Therefore, (m) strongly increases with A p . Relatively small reductions 
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us of rainfall at the transition between these two regimes may entail a dramatic increase in 

139 long-term soil salinization. Figure 2 also shows the threshold of soil salinity below which 

i« vegetation is practically unaffected (e.g., (C) < 2 dS/m) and the thresholds above which 

i4i regular (e.g. non-halophytic) vegetation is damaged (e.g., (C) > 2 dS/m). For coastal 

«2 areas soil salinization may occur even in relatively more humid regions, especially when 

143 rainfall events are not very intense. On the contrary, in continental regions only arid 

144 climates may begin to develop soil salinization (in the absence of irrigation and ground- 

145 water input). Indeed, through our model one can evaluate the risk of soil salinization in 
M 6 rain- fed agriculture just by estimating the typical salt inputs, total rainfall per year and 
147 the rainfall frequency. For example, a rain-fed crop in a semi-arid climate (e.g., rainfall 
us depth of 70 cm/yr) in a continental region risks salinization only when rainfall events are 

149 not very intense (e.g., 7" 1 < 0.4 cm or A p > 0.48 d _1 ). If the same crop is located in a 

150 coastal area, salinization occurs for a wider range of rainfall parameters (e.g., 7" 1 < 1 cm 

151 or A p > 0.18 d" 1 ). 

152 The solution (7) may be used in conjunction with soil moisture statistics to obtain a full 

153 characterization of the salt concentration in the root zone. Because one may safely assume 

154 that equations (1) and (3) are decoupled over short time scales, the soil moisture s(t) 

155 and the salt mass m(t) may be treated as statistically independent random variables. By 

156 observing that the salt concentration in the root zone is equal to C(t) = m{t)/nZ r s{t) and 

157 assuming s w ~ 0, we find the stationary probability distribution of the salt concentration 
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us p(C) as the quotient distribution of two independent random variables [Curtiss, 1941], 



The comparison between analytical solutions and numerical simulations (Figure 3) shows 

wo that the analytical solution reproduces reasonably well the pdf of the complete model, 

lei By integrating equation (8) from a given concentration value C* to infinity, one obtains 

162 the cumulative pdf of C, P(C*), which is the probability of having a salt concentration 

163 greater than a certain critical concentration value, C*, as a function of the soil-plant- 

164 atmosphere parameters. The inset of Figure 3 confirms the impact that climate change 

165 may have on soil salinity. Note, in particular, that such an impact is marked only for 
me semi-arid or drier climates (see Figure 2). For example with a reduction from A p = 0.2 

167 to A p = 0.15 d -1 , the probability of crossing C* = 6 dS/m is more than tripled. When 

168 coupled to a crossing analysis of concentration levels, the previous results may be used to 

169 evaluate the risk of plant salt stress. The analytical form of the results makes it suitable for 

170 computations of salinity risk at the global scale as a function of few measurable parameters, 

171 and facilitates their coupling with other models of long-term soil-plant biogeochemistry. 

4. Conclusions 

172 In this Letter we have presented an analytical approach to stochastic modeling of soil 

173 salinity, where the complexity of the problem is reduced by employing simplifying assump- 

174 tions that permit us to describe high-dimensional, unpredictable components via suitable 

175 random terms. By assuming time-averaged inputs of salt and instantaneous percolation 



P(C) 
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176 processes, a decoupling from soil moisture equation results in a simplified stochastic mass 

177 balance equation for the soil salt mass amenable to exact solution. 

178 Soil salinity statistics are obtained as a function of climate, soil and vegetation parame- 

179 ters. These can be combined with soil moisture statistics to obtain a full characterization 
wo of soil salt concentrations and the ensuing risk of primary salinization. 

181 This modeling framework can be extended to investigate additional salt inputs from 

182 irrigation and groundwater by modifying accordingly the average salt input parameter T 

183 and calculating the corresponding soil moisture pdfs (e.g. see Vervoort and Van Der Zee 

184 [2008] for groundwater inputs and Vico and Porporato [2009] for irrigation). 
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Figure 1. Comparison of soil moisture and salinity models: (a) Temporal evolution (equation 
(1)) of s(t), forced by intermittent rainfall (A p = 0.1 d" 1 and l/j P = 1.79 cm). The blue dashed 
line refers to the minimalist model, while the continuous red line is the complete numerical model 
(see text for details), (b) Temporal evolution of root-zone salt mass for the complete numerical 
model (red line) and the minimalist model (blue dashed line), (c) Temporal evolution of the 
corresponding specific salt concentration C(t) = m(t)/nZ r s(t) in the root zone for the same 
two cases of lb. We transform the unit of measure of C from mg/ (cm m 2 ) to dS/m, by using 
mg/ (cm m 2 ) = 10 _1 mg/1. The soil and vegetation parameters employed for the simulation of 
the complete model are those typical for a sandy-loam soil, while the free parameters of the 
minimalist model are si = 0.8, b = 0.6. In particular for both models we used n = 0.45, Z r = 30 
cm, s w = 0.1, ET max = 0.35 cm/d, Cr ~ 3 mg l" 1 and Aid = 54 mg d _1 m -2 (coastal area). 
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1- 2dS/m 

Continental region 
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Figure 2. Contour plot of the asymptotic mean concentration of salt (C) from the exact 
solution of (m) as a function of yearly rainfall depth and frequency. The values reported in 
the legend refer to the corresponding salt concentration values with respect to the average soil 
moisture (s) (for its analytical expression see Porporato et al. [2004]). The contour lines represent 
significant soil salinity values (1,2,..., 8 dS/m). The parameter \i has been calculated through 
equation (5); the others are as in Figure 1 for the coastal region, while for continental areas 
T « 6 mg d~ x mT 2 . 
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C [dS/m] 

Figure 3. Comparison between the analytical form of p(C) for the minimalist model obtained 
via equation (8) in which the free parameters (b and Si) have been fitted (solid line) and the 
numerical simulations of the corresponding complete model (dashed line). The soil and hydro- 
climatic parameters are as in Figure 1. In the inset: probability of exceeding a soil salinization 
critical threshold C* (Jq?° p(C)dC) as a function of the rainfall frequency A p . 
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